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INTRODUCTION 


This  annual  sununary  report  describes  research  performed  from  14  August  1992 
through  31  August  1993  with  support  from  ONR  under  grant  N00014-89-J-1003. 
Three  main  projects  are  discussed  in  this  report: 

I.  Nonlinear  Rayleigh  Waves 

II.  Pulsed  Finite  Amplitude  Sound  Beams 

III.  Finite  Amplitude  Sotmd  in  Waveguides 

Contributions  to  these  projects  were  made  by  the  following  individuab: 

Senior  Personnel 

•  M.  F.  Hamilton,  principal  investigator 

•  Yu.  A.  n’insky,  visiting  scientist 

•  E.  A.  Zabolotskaya,  visiting  scientist 
Graduate  Students 

•  M.  A.  Averkiou,  Ph.D.  student  in  M&chanical  Engineering 

•  E.  Yu.  Knight,  M.A.  student  in  Physics 

•  Y.-S.  Lee,  Ph.D.  student  in  Mechanical  Engineering 

•  G.  D.  Meegan,  Ph.D.  student  in  Physics 

•  D.  J.  Shull,  M.S.  student  in  Electrical  Engineering 

•  T.  W.  VanDoren,  Ph.D.  student  in  Medianical  Engineering 

Although  Professor  Il’insky  was  in  Moscow  during  the  entire  period  covered  by 
this  report,  he  was  instrumental  in  the  research  performed  on  nonlinear  Rayleigh 
waves  as  a  result  of  his  visit  to  Aristin  during  the  1991-92  academic  year.  Professor 
Il’insky  was  consulted  frequently  by  our  group  during  the  past  academic  year. 
Dr.  Zabolotskaya  is  on  leave  from  the  General  Physics  Institute  in  Moscow. 

The  main  source  of  firumdal  support,  in  addition  to  that  provided  by  ONR, 
has  been  the  David  and  Lucile  Packard  Foundation  Fellowship  for  Science  and  En¬ 
gineering.  Computing  resources  were  provided  by  The  University  of  Texas  System 
Center  for  High  Performance  C' — nuting. 

The  following  maniiscripts  1  abstracts,  which  describe  work  supported  at 
least  in  part  by  ONR,  have  beei  .mblished  (or  submitted  for  publication)  since 
14  August  1992. 
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Refereed  Journals 


•  Yu.  A.  n'inskii  and  E.  A.  Zabolotskaya,  “Cooperative  radiation  and  scattering 
of  acoustic  waves  by  gas  bubbles  in  liquids,”  J.  Acoust.  Soc.  Am.  92,  2837- 
2841  (1992). 

•  M.  F.  Hamilton,  “IVansient  axial  solution  for  the  reflection  of  a  spherical 
wave  from  a  concave  ellipsoidal  mirror,”  J.  Acoust.  Soc.  Am.  93,  1256-1266 
(1993). 

•  C.  M.  Darvennes  and  M.  F.  Hamilton,  “Additional  remarks  on  parametric 
reception  near  a  reflecting  plane,”  J.  Acoust.  Soc.  Am.  93,  3507-3510  (1993). 

•  M.  F.  Hanxilton,  Yu.  A.  Il’insky,  and  E.  A.  Zabolotskaya,  “On  the  existence 
of  stationary  nonlinear  Rayleigh  waves,”  J.  Acoust.  Soc.  Am.  93,  3089-3095 
(1993). 

•  D.  J.  Shull,  M.  F.  Hamilton,  Yu.  A.  Il’insky,  and  E.  A.  Zabolotskaya,  “Har¬ 
monic  generation  in  plane  and  cylindrical  nonlinear  Rayleigh  waves,”  J.  Acoust 
Soc.  Am.  94,  418-427  (1993). 

•  M.  A.  Averkiou,  Y.-S.  Lee,  and  M.  F.  Hamilton,  “Self-demodulation  of  am¬ 
plitude  and  frequency  modulated  pulses  in  a  thermoviscous  fluid,”  J.  Acoust. 
Soc.  Am.  (in  press). 

Conference  Proceedings 

s  M.  F.  Hamilton,  “A  transient  solution  for  the  axial  pressure  field  of  a  spark- 
source  lithotripter,”  Proceedings  of  the  14tb  International  Congress  on  Acous¬ 
tics,  edited  by  P.  Li  (Beijing,  China,  1992),  Vol.  1,  Paper  A5-1. 

•  M.  A.  Averkiou,  Y.-S.  Lee,  and  M.  F.  Hamilton,  “Self-demodulation  revis¬ 
ited,”  Advances  in  Nonlinear  Acoustics,  edited  by  H.  Hobsek  (World  Scien¬ 
tific,  Singapore,  1993),  pp.  251-256. 

•  D.  J.  Shull,  M.  F.  Hamilton,  and  E.  A.  Zabolotskaya,  “Nonlinear  Rayleigh 
wave  beams,”  Advances  in  Nonlinear  Acoustics,  edited  by  H.  Hobsek  (World 
Scientific,  Singapore,  1993),  pp.  496-501. 

Oral  Presentation  Abstracts 

•  D.  J.  Shull,  M.  F.  Hamilton,  Yu.  A.  Il’insky,  and  E.  A.  Zabolotskaya,  “Har¬ 
monic  interactions  in  plane  and  cylindrical  nonlinear  Rayleigh  waves,”  J. 
Acoust.  Soc.  Am.  92,  2358  (1992). 

•  M.  F.  Hamilton,  Yu.  A.  Il’insky,  and  E.  A.  Zabolotskaya,  “On  the  existence  oi 
stationary  nonlinear  Rayleigh  waves,”  J.  Acoust.  Soc.  Am.  92,  2358  (1992). 
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•  M.  F.  Hamilton,  Yu.  A.  Il’insky,  and  E.  A.  Zabolotskaya,  “Rayleigh  wave 
nonlinearity,”  J.  Acoust.  Soc.  Am.  93,  2384  (1993). 

These* 

•  C.  E.  Bruch,  “Second  harmonic  generation  in  Pekeris  waveguides,”  M.A. 
Thesis,  The  University  of  Texas  at  Austin  (December  1992). 

•  D.  J.  Shull,  “Harmonic  generation  in  plane,  cylindrical,  and  diffracting  non¬ 
linear  Rayleigh  waves,”  M.S.  Thesis,  The  University  of  Texas  at  Austin  (May 
1993). 

•  T.  W.  VanOoren,  “Propagation  of  finite  amplitude  sound  in  multiple  wave¬ 
guide  modes,”  Ph.D.  Dissertation,  The  University  of  Texas  at  Austin  (August 
1993). 

I.  Nonlinear  Rayleigh  Waves 

Portions  of  the  work  described  in  this  section  were  performed  by  Yu.  A.  D’insky, 
E.  Yu.  Knight,  G.  D.  Meegan,  D.  J.  Shull,  and  E.  A.  Zabolotskaya.  Three  projects 
are  described  below:  diffracting  nonlinear  Rayleigh  wave  beams,  pulsed  nonlinear 
Rayleigh  waves,  and  local  and  nonlocal  nonlinearity  in  Rayleigh  waves.  All  three 
projects  are  based  on  the  theoretical  model  developed  by  Zabolotskaya^  for  the 
propagation  of  nonlinear  Rayleigh  waves  in  isotropic  solids. 

A.  Diff'racting  Nonlinear  Rayleigh  Wave  Beams 

This  project  is  completed,  and  it  was  reported  first  in  Shull’s  masters  thesis.^ 
Shull  graduated  with  a  masters  degree  in  Electrical  Enpneering  in  May  1993.  The 
following  disctission  is  mainly  excerpted  from  the  proceedings’  that  accompanied 
the  presentation  of  this  material  at  the  13th  International  Symposium  on  Nonlinear 
Acoustics  (Bergen,  Norway,  Jvme-July  1993).  An  article  is  currently  in  preparation 
for  submission  to  Journal  of  the  Acoustical  Society  of  America. 

Linear  theory  for  diffracting,  monofrequency,  Rayleigh  wave  beams  is  well 
established.^  Here,  nonlinear  effects  that  produce  harmonic  generation  and  wave¬ 
form  distortion  in  diffracting  Rayld^  wave  beams  are  taken  into  account.  The 
isotropic  elastic  half  space  that  supports  the  Rayleigh  wave  occupies  the  region 
z  <  0.  The  source  is  in  the  plane  x  =  0  and  it  radiates  in  the  -i-x  direction.  Of 
interest  here  is  the  effect  of  finite  source  width  in  the  y  direction,  which  leads  to 
diffraction  in  the  (x,  y)  plane. 

For  simplicity,  all  equations  presented  below  are  in  dimensionless  form.  A 
complete  list  of  definitions  for  all  parameters  needed  to  perform  the  calculations 
may  be  found  elsewhere.’*’  Only  parameters  that  are  necessary  for  interpretation 
of  the  resxilts  are  provided. 
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The  analysis  applies  to  the  honzontal  (v^)  and  vertical  (v^)  components  of  the 
particle  velocity  vector.  The  components  are  written  in  the  dimensionless  forms 
14  =  Vgfvog  and  V,  =  u,/uo„  where  uq,  and  uq,  are  characteristic  source  velocity 
amplitudes,  and  the  following  Fourier  series  expansions  are  introduced: 

=  v4x,Y)u,n(Z)e-'^^  ^*  =  5  E  v;(x,r)c;,n(Z)c-‘^  (i) 

where  X  =  x/xq  is  an  axial  coordinate  in  terms  of  the  Rayleigh  distance  zq  = 
hid?  12,  ko  —  uq/cr  is  a  wavenumber  in  terms  of  the  fimdamental  angular  frequency 
uo  of  the  sotirce  and  the  small  signal  Rayleigh  wave  speed  ca,  a  is  a  characteristic 
source  dimension,  Y  =  y/a  is  a  dimensionless  transverse  coordinate,  Z  =  koz  is 
a  dimensionless  depth  coordinate,  and  r  =  u;o(t  —  x/cr)  is  a  retarded  time.  The 
functions  Usn  &nd  Utn  describe  the  depth  dependencies,  according  to  linear  theory, 
for  the  individtial  harmonic  components.  The  following  set  of  coupled  equations 
can  then  be  derived,  within  the  parabolic  approximation,  for  the  finite  amplitude 
behavior  of  the  spectral  components  in  a  diffracting  beam:^ 

dV  1 

w  ^  ay? + =  -fl"’  E  (2) 

fUS^OO 

where  l  =  n  —  m,  =  ttnXo  is  an  absorption  parameter  in  terms  of  the  small  signal 
attenuation  coefficient  On  at  frequency  nte>o«  B  =  zo/7  is  a  dimensionless  source 
amplitude  expressed  in  terms  of  the  distance  7  that  characterizes  the  nominal  shock 
formation  distance  for  a  plane  wave,  and  R„^i  is  a  matrix  of  coupling  coefficients  for 
the  harmonic  interactions.  A  model  similar  to  Eq.  (2),  but  without  an  absorption 
term,  was  proposed  by  Parker.^ 

Asymptotic  quasilinear  solutions  of  Eq.  (2)  for  the  farfield  of  the  fundamen¬ 
tal  and  second  harmonic  components  have  been  derived  for  arbitrary  amplitude 
distributions  at  the  source.  The  results  are  similar  (apart  from  decay  rates)  to 
those  obtained  for  sound  beams  in  fluids,  and  shall  not  be  presented  here.  Instead, 
attention  is  focused  on  numerical  solutions  of  the  fully  nonlinear  system. 

To  integrate  Eq.  (2)  numerically  with  the  diffraction  term  included,  it  is  con¬ 
venient  to  introduce  a  coordinate  system  which  follows  the  divergent  spreading  of 
the  beam.*  The  appropriate  transformation  for  a  cylindrically  spreading  beam  is 
defined  as  follows:*** 

r,(jr,F)  =  v.(jr,y)v^rMexp(-|i^) 

where  ^  is  a  parameter  that  controls  the  divergence  (S  >  0)  or  convergence  (S  <  0) 
of  the  original  coordinate  system.  Elimination  of  Vn  and  Y  in  Eq.  (2),  in  favor  of 
Vn  and  y,  yields 

W  +  4.n(X+<)»^’  +  "  vTTI  (4) 
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The  resiilts  presented  below  were  obtained  by  solving  Eq.  (4)  with  finite  difference 
methods  (following  minor  modifications)  that  are  commonly  used  to  integrate  the 
spectral  form  of  the  KZK  equation.*® 

The  results  were  calculated  for  a  monofrequency  source  with  a  uniform  rectan¬ 
gular  amplitude  distribution  as  a  function  of  the  transverse  coordinate  V ; 

v;(o,y)  =  i  |r|<i  (5) 

=  0  |r|  >  1 

Vn(0,r)  =  0  n>l 

This  sotirce  condition  is  a  common  model  for  interdigital  transducers  that  are  used 
on  surface  acoustic  wave  devices.^  Investigations  based  on  other  source  conditions, 
including  foctised  and  asymmetric  sources,  are  reported  elsewhere.^*®  The  values 
Ai  =  10~^  and  B  =  I  were  used  in  all  computations,  and  the  attenuation  co¬ 
efficients  depend  quadratically  on  frequency  (^4^  3=  All  niimerical  results 

presented  below  are  for  the  free  stirface  Z  =  0,  and  they  are  based  on  a  coefficient 
matrix  Rmi  that  was  evaluated  with  parameters  for  steel.**’ 

Presented  in  Fig.  1  are  (a)  axial  propagation  ctirves  and  (b)  beam  patterns  for 
the  lowest  four  harmonic  components  in  a  Rayleigh  wave  beam  of  finite  amplitude. 
The  results  are  qualitatively  similar  to  those  for  finite  amplitude  sotmd  that  is 
radiated  into  a  viscous  fluid  by  a  uniform  circular  source.*  In  particular,  efficient 
harmonic  generation  occurs  mainly  in  the  nearfield  (X  <  1),  and  the  beam  pattern 
for  the  nth  harmonic  component  possesses  n  times  as  many  sidelobes  as  are  pre¬ 
dicted  by  linear  theory.  The  beam  pattern  for  the  fundamental  component  has  a 
flatter  main  lobe  and  higher  sidelobes  than  are  predicted  by  linear  theory  (because 
of  shock  formation  in  the  main  lobe),  which  is  again  consistent  with  the  results  for 
finite  amplitude  sotmd  in  fluids. 

Axial  waveforms  for  the  horizontal  velocity  component,  at  three  distances  from 
the  source,  are  presented  in  Fig.  2.  At  X  =  1.5,  the  sharpening  of  the  peak  and  the 


Figure  1:  (a)  Axial  propagation  curves  (y  =  0)  and  (b)  beam  patterns  at  X  =  10 
for  the  lowest  four  harmonic  components  in  a  Rayleigh  wave  beam. 
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rounding  of  the  trough  are  characteristic  of  the  effect  of  diffraction  on  waveforms 
in  fluids.  However,  the  waveform  asymmetry  in  a  cylindrically  spreading  (two 
dimensional)  beam  is  less  than  in  a  spherically  spreading  (three  dimensional)  beam. 
The  waveform  {^symmetry  in  Fig.  2  is  even  less  noticeable  at  =  3,  where  a  cusped 
sawtooth  profile,  similar  to  that  predicted  for  plane  waves,  has  formed.  Note  the 
approximately  7r/4  phase  advance  in  the  positions  of  the  shocks,  which  agrees  with 
linear  diffraction  theory  for  the  phase  of  the  fundamental  component  in  the  farfield. 


Finally,  transverse  amplitude  distributions  for  the  lowest  three  harmonic  com¬ 
ponents  in  the  nearfield  of  the  beam,  at  =  0.18  (which  is  the  location  a£  the 
last  axial  minimtim  in  the  neaifield  of  the  fimdamental  component),  are  presented 
in  Fig.  3.  The  distribution  for  the  fimdamental  component  is  very  close  to  that 
predicted  by  linear  theory.^ 


Figure  3:  Transverse  harmonic  amplitude  distributions  in  the  nearfield  of  a  beam 
(X  —  0.18).  A  different  vertical  scale  is  used  for  each  of  the  harmonics. 


B.  Local  and  Nonlocal  Nonlinearity  in  Rayleigh  Waves 

Results  from  this  work  were  presented^^  at  the  125th  Meeting  of  the  Acoustical 
Society  of  America  in  Ottawa,  Canada  (May  1993).  An  article  is  currently  in 
preparation  for  submission  to  Journal  of  the  Acoustical  Society  of  America. 

Previous  theoretical  investigations  of  waveform  distortion  and  shock  forma¬ 
tion  in  finite  amplitude  Rayleigh  waves  are  based  on  periodic  waveforms,  mainly 
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waveforms  that  are  siniisoidal  at  the  source.  Theoretical  investigations  of  pulsed 
nonlinear  Rayleigh  waves,  e.g.,  those  which  are  exploited  in  nonlinear  signal  pro¬ 
cessing  applications  of  surface  acoustic  wave  (SAW)  devices,  are  based  on  simple 
quasilinear  modeb  of  the  form^^ 

i,3(t)  =  f,  t  -  dUdlT  (6) 

Here  Vx(t)  and  t^(t)  are  signals  whose  weak  nonlinear  interaction  produces  a  sec¬ 
ondary  signal  V3(f)  at  a  given  location,  h(ti,t3)  is  an  impulse  response,  and  T  b 
an  interval  that  includes  the  duration  of  the  interaction  (or,  alternatively,  it  b  an 
integer  number  of  periods  in  a  time  harmonic  system).  In  modeb  such  as  Eq.  (6), 
h(ti,t2)  b  simply  a  delta  function  that  accoimts  for  time  delays  associated  with 
the  interacting  waves. 

Several  limitations  of  Eq.  (6)  can  be  noted.  First,  the  theoretical  model  b 
quasilinear  and  can  thus  accoimt  only  for  weakly  nonlinear  effects,  such  as  stun 
and  difference  frequency  generation  in  harmonic  waves,  or  products  which  can  be 
associated  with  convolutions  or  correlations  of  the  pulsed  signab  vj  and  V2.  It 
cannot  account  for  the  nonlinear  distortion  of  a  Rayleigh  as  it  propagates,  and  it 
cannot  account  for  shock  formation.  Second,  the  impulse  response  h(<i,  tj)  b  based 
on  empirical  results — no  analytical  expressions  for  h(ti,t2)  have  appeared  in  the 
literature.  Third,  the  assiunption  that  h(t|,t2)  is  proportional  to  delta  functions 
is  equivalent  to  the  assumption  that  the  nonlinearity  in  Rayleigh  waves  b  entirely 
local.  However,  it  is  known  that  nonlinearity  in  Rayleigh  waves  b  nonlocal,  as 
pointed  out  first  by  Parker  and  Talbot.*^ 

Zabolotskaya^  showed  that  the  evolution  equation  for  the  horizontal  velocity 
component  at  the  siuface  an  isotropic  solid  may  be  written  in  the  form 

where  v(x,r)  is  the  horizontal  velocity  component,  z  is  dbtance,  r  ^  t  —  xfcR  is 
a  retarded  time,  L  is  the  kemal  of  the  integral  operator,  and  C  is  a  dimensionless 
constant.  In  her  original  paper,^  the  kernel  L  was  expressed  as  the  Fourier  series 

i(‘l .  ‘j)  =  E  !"•  +  'I  A.!*''””"  (8) 

mj 

where  is  the  same  matrix  that  appears  in  Eq.  (2). 

During  the  past  year,  it  was  shown  that  L  can  be  expressed  in  the  explicit 
time-domain  form 

where  A  is  a  dimensionless  constant,  B  is  dimensionless  function  of  the  ratio  U/t^, 
and  6  is  the  Dirac  delta  fxmction.  Explicit  expressions  for  A  and  B  (as  well  as 
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C),  in  terms  of  second  and  third  order  elastic  constants  of  the  material,  have  been 
derived.  Substitution  of  Eq.  (9)  in  Eq.  (7)  yields 


^  =  £  \a—  4.  a  //  p  /  r  -  \  t;(x,  r^)  dr' dr'' 

dx  cjf  dr  dr  JJt  \r  —  (r  —  +  (r  —  r")^ 


(10) 


Setting  J?  =  0  in  Eq.  (10)  produces  an  equation  identical  in  form  to  that  for  plane 
soimd  waves  of  finite  amplitude  in  a  lossless  fiuid.  Nonlinearity  in  fiuids  is  entirely 
local:  finite  amplitude  distortion  at  an  instant  on  the  waveform  is  determined  (in 
the  preshock  region)  by  the  value  of  dv^/dr  at  that  same  instant.  The  first  term 
on  the  right-hand  side  of  Eq.  (9)  is  thus  associated  with  local  nonlinear  effects.  In 
contrast,  nonlinearity  due  to  the  integral  in  Eq.  (10),  at  a  given  instant  r,  depends 
on  the  value  of  v  at  all  instants  r.  The  second  term  on  the  right-hand  side  of  Ek|.  (9) 
is  thus  associated  with  nonlocal  nonlinear  effects.  The  nonlocal  nonlinearity  in 
Rayleigh  waves  is  what  distinguishes  their  finite  amplitude  properties  from  those 
of  sotmd  waves  in  fluids. 

Comparison  of  Eqs.  (6)  and  (10)  illustrates  the  limitations  of  Eq.  (6)  that  were 
discussed  above.  Equation  (10)  describes  the  evolution  of  the  signal  throughout  its 
path  of  propagation,  not  just  at  a  single  point;  it  accotmts  for  waveform  distortion 
and  shock  formation;  it  is  derived  rigorously  in  terms  cf  second  and  third  order 
elastic  constants  of  the  material;  and  it  accounts  for  nonlocal  as  well  as  local 
nonlinear  effects. 

Althou^  Eq.  (10)  is  an  explicit  evolution  equation  for  nonlinear  Rayleigh  waves, 
and  it  provides  substantial  insight  into  the  natiire  of  the  nonlinearity,  its  form 
is  not  the  most  convenient  for  performing  munerical  calculations  (e.g.,  such  as 
those  required  to  describe  waveform  distortion,  as  depicted  in  Fig.  2).  Instead,  a 
differential  evolution  evolution  which  is  currently  imder  development  [Eq.  (4)  of  the 
previous  annual  summary  report^^]  provides  a  more  practical  means  of  performing 
time  domain  calculations  for  nonlinear  Rayleigh  waves. 


C.  Pulsed  Nonlinear  Rayleigh  Waves 

The  following  results  for  ptilsed  Rayleigh  waves  were  obtained  by  Knight,  who 
entered  the  M.A.  program  in  the  Phyacs  Department  in  fall  1992.  Knight  receives 
her  salary  support  from  Los  Alamos  National  Laboratory. 

As  noted  above,  there  are  no  reported  theoretical  predictions  of  nonlinear  distor¬ 
tion  and  shock  formation  in  pulsed  Rayleigh  waves.  In  fact,  munerical  simulations 
of  even  periodic  Rayleigh  waveforms,  well  beyond  the  point  of  shock  formation,  have 
been  calculated  only  on  the  basis  of  Zabolotskaya's  model  equations.  There  are 
presently  three  approaches  that  may  be  followed  to  investigate  pulsed  nonlinear 
Rayleigh  waves  ntunerically.  One  is  to  construct  repeated  sequences  of  pulses  with 
harmonic  series  and  then  solve  an  equation  similar  to  Eq.  (2).  Equation  (7)  can 
also  be  used,  but  the  convolution  is  very  time  consuming.  Finally,  the  differential 
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evolution  equation  which  is  discussed  in  the  previous  annual  summary  report*^  can 
be  used  (this  work  is  currently  in  progress).  Here,  the  first  approach  is  used  to 
obtain  numerical  results  for  plane  wave  pulses.  Equation  (2)  is  rewritten  in  terms 
of  different  dimensionless  quantities: 

dV  °° 

+  AnVn  =  (11) 

mss^oo 

where  X  =  xfx  and  An  =  oinX.  The  distance  and  attenuation  parsuneters  are  now 
scaled  according  to  7,  the  nominal  shock  formation  distance  in  steel  at  the  center 
frequency  of  the  pulse.  The  center  frequency  is  taken  to  be  the  iVth  harmonic  in 
the  Fourier  series  expansion.  Waveforms  are  again  obtained  from  Eqs.  (1).  Care 
must  be  taken  that  the  repetition  rate  of  the  pulses  described  by  the  Fourier  series 
is  sufficiently  low  that  successive  pulses  do  not  interact  with  each  other. 

Shown  in  Fig.  4  are  predictions  for  the  propagation  of  a  short  tone  burst  (center 
frequency  u^o)  with  a  Gaussian  envelope  function,  for  An  =  0.01.  The  center 
frequency  of  the  source  waveform,  u;o,  is  iV  =  20  times  the  fundamental  frequency 
in  the  Fourier  series  representation.  Replicas  of  the  pulses  shown  in  Fig.  4  are 

repeated  at  u;oT/2ir  =  20n,  where  n  =  ±1,±2, _  The  number  of  harmonics 

retained  in  the  computations  was  500,  which  corresponds  to  25  harmonics  of  the 
center  frequency  at  the  source. 

The  horizontal  and  vertical  velocity  waveforms  are  shown  (as  functions  of  dis¬ 
tance  X)  in  the  second  and  third  columns  of  Fig.  4,  respectively.  In  the  first  col- 
unm,  for  comparison,  are  predictions  for  the  propagation  of  sotmd  in  a  fluid  (based 
on  the  Burgers  equation),  at  distances  where  the  distortion  is  similar  to  that  in 
the  horizontal  component  of  the  Rayleigh  wave  particle  velocity.  The  ranges  at 
which  the  calculations  for  fliiids  are  compared  with  the  calculations  for  Rayleigh 
waves  are  based  on  comparisons  of  the  shock  formation  distances.  Speciflcally, 
waveforms  are  compared  at  distances  that  are  the  same  when  normalized  by  the 
distance  where  shock  formation  first  occurs  in  an  initially  sinusoidal  wave  at  the 
center  frequency  of  the  sovirce.  The  dimensionless  absorption  coefficients  are  also 
adjusted  accordingly. 

Comparison  of  the  waveforms  in  the  first  two  columns  reveals  that  the  propaga¬ 
tion  a  Rayleigh  wave  pulse  resembles  somewhat  the  propagation  of  an  acoustic 
pulse  in  a  fluid.  One  main  distinction,  as  in  the  case  of  periodic  waves, is  the 
cusping  near  shocks  in  the  Rayleigh  waves.  In  contrast,  sound  waves  in  fltiids 
have  more  of  a  sawtooth  profile  (e.g.,  at  X  =  10).  Beyond  the  shock  wave  region, 
both  waveforms  experience  self-demodulation  and  tend  toward  a  wave  shape  that 
is  more  closely  related  to  the  envelope  of  the  source  waveform  than  to  the  center 
frequency.  Comparison  of  the  waveforms  at  X  =  400  and  X  =  1200  reveals  that 
the  absorption  in  this  example  is  sufficiently  weak  that  the  demodulated  waveforms 
are  begiiming  to  form  shocks.  Moreover,  the  ‘‘N  wave”  predicted  in  the  fluid  at 
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Figure  4:  Horizont^d  velocity  waveforms  (second  column),  vertical  velocity  wave¬ 
forms  (third  column),  and  corresponding  frequency  spectra  (fourth  column)  for  a 
nonlinear  Rayleigh  wave  pulse,  as  functions  of  distance.  In  the  first  column  are 
comparable  results  for  sotmd  waves  in  fluids.. 
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X  —  1200  has  twice  the  amplitude  of  the  horizontal  component  of  the  Rayleigh 
wave,  and  the  head  and  tail  shocks  have  already  begun  to  move  away  from  each 
other. 


01234S80  1  2  3 


(O/COo  (O/ttg 

Figure  5:  Comparison  of  frequency  spectra  for  pulses  similar  to  those  in  Fig.  4,  but 
with  higher  attenuation.. 

In  Fig.  5  are  shown  comparisons  of  frequency  spectra  for  Rayleigh  waves  (solid 
lines)  and  for  sound  waves  (dashed  lines)  in  fluids,  again  at  distances  that  are 
matched  according  to  the  scheme  described  above.  Here,  the  absorption  is  ten 
times  greater  than  in  Fig.  4  {As  =  0.1).  At  X  =  2,  the  generation  of  higher 
harmonics  in  Rayleigh  waves  and  sound  waves  is  very  similar,  both  in  terms  of 
spectral  amplitudes  and  bandwidths.  Note,  however,  that  the  generation  of  ‘‘dif¬ 
ference  frequencies”  (i.e.,  spectral  components  below  the  frequency  band  of  the 
primary  wave)  is  more  efficient  in  sound  waves  than  in  Rayleigh  waves.  [This  ef¬ 
fect  can  be  demonstrated  analytically  by  comparing  quasilinear  solutions,  for  sotmd 
waves  and  Rayleigh  waves,  for  the  case  of  s\un  and  difference  frequency  generation.] 
The  spectra  in  Fig.  5  thus  assist  interpretation  of  the  waveforms  in  the  flrst  two 
columns  of  Fig.  4.  Up  to  A*  =  10  in  Fig.  4,  the  processes  of  distortion  and  shock 
formation,  which  are  determined  mainly  by  higher  harmonic  generation,  are  very 
similar  in  sotmd  waves  and  Rayleigh  waves,  ki  X  =  40,  the  stronger  difference 
frequency  generation  in  the  fluid  can  be  noticed.  Comparison  of  the  sound  wave 
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at  =  200  and  the  Rayleigh  wave  9X  X  —  400  reveals  that  the  sound  wave  in  this 
example  demodulates  in  approximately  half  the  distance  required  for  the  Rayleigh 
wave,  and  the  finsd  waveform  in  the  fluid  has  nearly  twice  the  amplitude  as  the 
Rayleigh  wave. 

D.  Nonlinear  Stoneley  waves 

This  research  is  being  performed  by  Meegan,  who  is  supported  by  ONR  through  the 
AASERT  program.  Meegan  entered  the  Ph.D.  program  in  the  Physics  Department 
in  faU  1992. 

The  main  focus  of  Meegan’s  research  is  the  development  of  a  theoretical  model 
for  nonlinear  Stoneley  waves.  The  goal  is  to  extend  the  model  equation  for  Rayleigh 
waves,  Eq.  (11),  to  account  for  the  spectral  interactions  in  Stoneley  waves.  In  prin¬ 
ciple,  the  method  of  accomplishing  this  task  is  straightforward.  The  Hamiltonian 
formalisr^  employed  by  Zabolotskaya^  to  derive  Eq.  (11)  for  Rayleigh  waves  is  eas¬ 
ily  adapted  for  application  to  Stoneley  waves.  The  main  difficulty  is  the  amount  of 
algebra  reqtiired  to  derive  the  matrix  Rni,  which  is  more  cumbersome  for  Stoneley 
waves  than  for  Rayleigh  waves  (even  the  algebra  for  Rayleigh  waves  is  exceedin^y 
tedious  to  accomplish  by  hand).  The  plan  is  to  use  the  Mathematica  symbolic 
computation  software  to  overcome  the  obstacle  presented  by  the  algebra.  The  first 
step,  the  development  a  Mathematica  program  which  reproduces  Zabolotskaya’s 
results  for  Rayleigh  waves, ^  is  completed.  The  relevant  material  parameters  re- 
qmred  to  analyze  nonlinear  Stoneley  waves  have  been  identified,  and  therefore  the 
remaining  task  is  to  modify  the  Mathematica  program  to  calculate  the  nonlinearity 
matrix. 

A  second  objective  of  Meegan’s  research  is  to  obtain  measurements  of  nonlin¬ 
ear  Rayleigh  waves  that  would  complement  recent  theoretical  predictions.^’’  Otir 
laboratory  was  not  equipped  previously  to  either  generate  or  measure  Rayleigh 
waves,  and  therefore  time  was  devoted  to  development  of  the  required  facilities. 
Pinducer  receivers  have  been  purchased,  and  piezoelectric  wedge  transducers  have 
been  obtained  for  tise  as  both  sources  and  receivers.  Preliminary  measurements  in 
aluminum  have  indicated  the  existence  of  harmonic  generation  due  to  finite  am¬ 
plitude  propagation.  However,  poor  signal-to-noise  ratios  are  encoimtered  even  at 
the  second  harmonic  frequencies,  and  therefore  more  efficient  source  and  receiver 
configurations  are  being  considered. 
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II.  Pulsed  Finite  Amplitude  Sound  Beams 

Portions  of  the  work  described  in  this  section  were  performed  by  M.  A.  Averkiou 
and  Y.-S.  Lee.  Lee  successfully  defended  his  Ph.D.  dissertation  in  the  Mechanical 
Engmeering  Department  in  July  1993.  His  dissertation  describes  the  development 
of  a  computer  code  for  solving  the  KZK  equation  for  pulsed  sound  beams  of  finite 
amplitude  in  thermoviscous  fiuids.  After  he  makes  minor  revisions  in  the  disserta¬ 
tion,  Lee  will  be  eligible  to  graduate  in  December  1993. 

A.  Uniformly  Valid  Axial  Solution  for  Self-Demodulation 

This  project,  performed  by  Averkiou  and  Lee,  is  completed  and  is  reported  in  two 
publications:  in  the  Proceedings  of  the  13th  International  Symposium  on  Nonlinear 
Acoustics^’  (Bergen,  Norway,  Jvme-July  1993);  and  in  an  article  that  is  pending 
publication  in  Journal  of  the  Acoustical  Society  of  America}^  The  following  dis¬ 
cussion  follows  the  former  publication.^^ 

Our  investigation  of  the  classic  problem  of  self-demodulation  was  described  in  a 
previous  annual  report  where  comparison  of  experiment  was  made  with  numerical 
predictions.^^  In  the  process  of  completing  that  work,  a  simple  yet  very  accurate 
analytic  solution  was  developed  for  the  field  along  the  axis  of  the  sound  beam.  The 
source  is  assumed  to  be  a  circular  piston  that  vibrates  with  uniform  amplitude  in 
the  plane  <7  —  0,  i.e., 


P(0,p,r)  = /(t)  p<1 

=  0  p>l 


(12) 


where  P{<r,  p,  r)  is  a  dimensionless  pressure,  a  is  dimensionless  range  in  terms  d 
the  Rayleigh  distance,  p  is  a  dimensionless  coordinate  transverse  to  the  axis  of  the 
beam,  normalized  by  the  piston  radiiis,  and  r  is  a  dimensionless  retarded  time. 
The  time  dependence  /(r)  is  defined  in  terms  of  an  amplitude  modulation  E{t) 
and  phase  modulation  ^r)  of  a  sinusoidal  oscillation,  as  follows: 

/(r )  =  E{t)  dn[r  +  ( 13) 


The  dimensionless  expression  for  the  instantaneous  angular  frequency  is  thus  D(r)  = 
1  +  d<f>/dT,  Now  assume  that  E(t)  and  ^(r)  are  slowly  varying  functions  in  com¬ 
parison  with  sinr.  Combining  various  results  of  others,^*' we  obtain  for  A  ^  1 
(absorption  length  less  than  the  Rayleigh  distance)  and  F  =  1  (where  F  is  the 

Gol’dberg  number),  the  following  solution  for  the  axial  waveform:^* 


P(<7,0,r)  = 


yJ^itAa 


where  the  asterisk  indicates  convolution  with  respect  to  r. 
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The  first  two  terms  within  the  square  brackets  of  Eq.  (14)  represent  the  cen¬ 
ter  wave  and  edge  wave,  respectively,  associated  with  the  axial  field  of  the  primary 
beam  radiated  by  a  \miform  piston  source.  The  convolution  of  these  two  terms  with 
the  exponential  dissipation  function  outside  the  brackets  accovints  for  thermovis- 
cous  attenuation  of  the  primary  beam,  which  was  investigated  recently  by  FV0ysa, 
N2ize  Tjotta,  and  Tj0tta.^^  The  third  term  in  the  square  brackets  is  the  farfield 
residt  for  the  low  frequency,  secondary  pressure  waveform  produced  by  the  distor¬ 
tion  of  the  primary  beam.  With  ^  =  const,  the  third  term  reduces  to  the  classical 
result  of  Berktay.**  The  factor  d^/dr  accounts  for  the  effect  of  phase  modulatitm, 
which  was  considered  by  Gurbatov,  Demin,  and  Malakhov.**  Convolution  of  the 
third  term  with  the  dissipation  function  approximates  the  effect  of  thermoviscous 
attenuation  on  the  secondary  wave,  in  accordance  with  the  approach  followed  by 
Cervenka  and  Alais.*° 

The  accuracy  of  Eq.  (14)  can  be  tested  agiunst  numerical  results  from  the  finite 
difference  solution  developed  by  Lee.*^  To  perform  this  comparison  we  consider  the 
limiting  values  >1  =  1  and  F  =  1  (agreement  is  expected  to  improve  for  larger  A 
or  smaller  F).  For  the  sotirce  waveform  we  let  E(r)  =  exp[— (r/25ir)®]  and  ^(r)  = 
— 5sin(r/25),  which  produces  a  pulse  with  sinusoidal  frequency  modulation.  The 
dimensionless  frequency  of  the  carrier  wave  thus  increases  from  Q  =  0.8  in  the 
center  of  the  pulse  to  D  1.2  at  either  end.  In  the  left  column  of  Fig.  6  is 
presented  a  comparison  of  time  waveforms  predicted  the  finite  difference  solution 
(solid  lines),  and  by  Eq.  (14)  (dashed  lines).  In  the  right  colunm  of  Fig.  6  are 
shown  the  corresponding  frequency  spectra  5,  normalized  to  yield  peak  amplitudes 
of  unity  at  the  soxirce  (<r  =  0).  The  waveform  and  firequency  spectrum  at  the 
sovirce  are  presented  in  the  top  row,  and  the  decibel  values  listed  in  each  frame 
indicate  the  level  relative  to  that  at  the  source.  Although  they  are  not  always  easily 
distinguished,  both  dashed  and  solid  lines  are  printed  in  each  frame  of  Fig.  6. 

Equation  (14),  despite  the  rather  stringent  conditions  under  which  it  was  de¬ 
rived,  is  in  good  agreement  with  the  numerical  solution  throughout  the  axial  field. 
At  <7  =  4,  the  main  effect  that  can  be  seen  is  due  to  thermoviscous  dissipation 
of  en«  I'gy  in  the  primary  beam,  which  causes  the  ends  of  the  pulse  to  attenuate 
more  rapidly  (because  o£  the  higher  instantaneous  frequency)  than  the  middle.  At 
<7  =  12,  the  demodulated  waveform  dominates  the  contributions  at  the  ends  (ff  the 
pulse,  but  the  primary  wave  still  dominates  the  middle.  The  final  demodulated 
waveform  is  achieved  at  <7  =  24.  The  positive  pressures  at  the  ends  of  the  pvilse 
are  due  primarily  to  amplitude  modulation,  and  the  negative  pressure  in  the  mid¬ 
dle  is  due  primarily  to  frequency  modulation.  Increasing  A  (with  respect  to  the 
value  of  unity  used  in  Fig.  6)  leads  to  predictions  by  Eq.  (14)  which  are  virtually 
indistinguishable  from  the  numerical  predictions.  Whereas  the  finite  difference  so¬ 
lution  reqtiires  computation  times  on  the  order  of  ten  minutes  on  a  super  computer, 
Eq.  (14)  requires  only  minutes  on  a  personal  computer. 
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Figure  6:  Left  column:  Comparison  of  the  axial  solution  (dashed  lines)  with  numer¬ 
ical  results  from  the  finite  difference  solution  (solid  lines)  for  a  frequency  modulated 
pulse  (<r  =  0),  with  A  =  F  =  1.  Right  column:  FVequency  spectre  of  the  waveforms 
in  the  left  coltuxm  (both  dashed  and  solid  lines). 


B.  Measurements  of  Pulsed  Sound  Beams  in  Water 

This  research  was  performed  by  Averkiou,  and  it  will  be  reported  at  the  fall  1993 
meeting  of  the  Acoustical  Society  of  America  in  Denver.  The  following  is  the 
abstract  that  was  submitted  for  that  meeting:^^ 

Measiuements  are  reported  for  finite  amplitude  acoustic  pulses  radiated 
by  plane  circular  pistons  in  water.  Pulses  with  center  frequencies  of  sev¬ 
eral  megahertz,  peak  sotmd  pressures  up  to  1  MPa,  durations  ranging 
from  approximately  two  to  twenty  cycles,  and  different  amplitude  and 
frequency  modulations  were  investigated.  Measurements  short  pulses 
were  made  very  near  the  sovurce,  where  the  center  wave  and  edge  wave 
can  be  separated.  P\ilse  envelope  distortion  that  accompanies  shock  for¬ 
mation  in  frequency  modulated  tone  bursts  is  demonstrated.  Acoustic 
saturation  of  pulsed  sotmd  beams  b  also  investigated.  All  measure¬ 
ments  are  compared  with  theoretical  predictions  obtained  from  a  com¬ 
puter  code  that  solves  the  KZK  equation  in  the  time  domain.^^  Very 
good  agreement  between  theory  and  experiment  is  obtained  both  in  the 
nearfield  and  the  farfield,  on  and  off  axis.  Artifacts  in  measurements  of 
waveforms  containing  shocks,  which  are  attributed  to  bandwidth  limi¬ 
tations  of  membrane  hydrophones,  axe  discussed. 

Representative  samples  of  the  results  are  shown  in  Fig.  7.  The  results  are  for 
radiation  from  a  1-inch  diameter  source  in  water  at  a  center  frequency  of  1  MHz. 
The  pulse  envelope  is  Gaussian,  and  tlw  equivalent  shock  formation  distance  for  a 
plane  wave  is  approximately  twice  the  Rayleigh  distance  of  the  source.  As  before, 
<7  =  z/zo  and  p  =  rfa.  The  first  three  rows  show  axial  waveforms  beginning  at 
the  source  (<r  =  0)  up  to  <7  =  2.9,  beyond  the  shock  fmmatiQn  distance.  The  last 
three  rows  show  wavefonns  across  the  beam  at  <r  =  2.9,  beginning  on  axis  (p  =  0) 
and  moving  out  to  p  5.  The  frequency  spectra  in  the  last  three  rows  reveal  the 
increase  in  directivity  as  a  function  of  increase  in  harmonic  number.  Agreement  of 
theory  and  experiment,  both  on  and  off  axis,  is  excellent. 
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Figure  7:  Waveforms  and  frequency  spectra  along  the  axis  (from  <7  =  0  to  =*  2.9) 
itnfl  transverse  to  the  axis  (from  p  =  0top  =  5)ofa  soimd  beam  in  water.  The 
solid  lines  are  calculated  and  the  dashed  lines  are  measured. 
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III.  Finite  Amplitude  Sound  in  Waveguides 

VanDoren  completed  his  research  on  this  project,  and  he  received  his  Ph.D.  in 
August  1993.  The  primary  source  of  VanDoren’s  salary  support  during  the  past 
two  years  was  a  Rockwell  Fellowship,  and  much  of  his  equipment  purchases  were 
fimded  by  the  Packard  Foundation.  Here  we  report  only  the  main  results  from  his 
dissertation. 

VanOoren  considered  the  analysis  of  a  waveguide  formed  by  a  horizontally  strat¬ 
ified  medium,  with  z  designating  the  axis  of  the  waveguide  and  z  the  transverse 
coordinate.  Acoustic  propagation  in  the  waveguide  is  modeled  by  the  following 
modified  form  dP  the  Westervelt  equation: 

_ _ MgX 

po4dt^  ^  ^ 

where  p  is  the  sound  pressure,  c(x)  accounts  for  weak  variation  of  the  sovmd  speed 
transverse  to  the  waveguide  axis,  and  the  right-hand  side  is  the  standard  nonlinear 
term  in  the  Westervelt  equation.  The  procedtire  for  solving  Eq.  (15)  is  an  extension 
of  the  method  outlined  in  Ref.  23.  It  is  assumed  that  the  solution  may  be  expressed 
in  the  form 

*  m  • 

where  are  eigenfunctions  that  describe  the  transverse  distributions  of  the  nor¬ 
mal  modes  (mode  m,  frequency  nw)  and  satisfy  boundary  conditions  at  the  walls 
(assxuned  parallel  to  the  z  axis),  pmn  account  for  amplitude  variations  along  the 
axis  of  the  waveguide,  and  are  axial  wavenumbers  obtained  from  linear  theory. 
If  the  wall  impedances  are  locally  reactive,  then  the  eigenfunctions  are  orthogo¬ 
nal,  and  the  following  apprcudznate  system  of  first-order  equations  can  be  obtained 
following  substitution  of  Eq.  (16)  into  Eq.  (15): 

Pm«(«)  +  a«nP-»n(*)  «  EE  E  +  E  (W 

’^miiPDVO  r  p  fml  gmt 

where  am*  are  attenuation  coefficients  that  were  introduced  ad  hoc,  the  fimction 
9mnptr  accou^ts  for  difference  frequency  generation,  and  accounts  for  sum 

frequency  generation.  The  summations  over  r  and  p  account  for  all  combinations 
of  modes  in  which  the  spectral  c<xnponents  can  interact,  and  the  sununation  over  g 
accounts  for  all  frequency  pairs  that  combine  to  yield  the  nth  harmonic  frequency 
component  on  the  left-hand  side.  Equation  (17)  was  solved  with  a  standard  fourth- 
order  Rimge-Kutta  method. 

A  main  advantage  of  the  above  solution  technique  is  that  the  parabolic  ap¬ 
proximation  is  avoided,  and  accurate  solutions  are  not  restricted  only  to  waves 


that  propagate  at  small  grazing  angles  (with  respect  to  the  axis  of  the  waveguide). 
In  addition,  the  modal  decomposition  that  is  inherent  in  the  solution  facilitates 
physical  interpretation  of  the  results.  An  alternative  and  currently  popular  model 
equation  for  investigating  finite  amplitude  sound  in  waveguides  is  the  Nonlinear 
Parabolic  E^uation^  (NPE).  However,  rather  than  make  comparisons  with  solu¬ 
tions  the  NPE,  we  compare  results  obtained  from  Eq.  (17)  with  finite  difference 
solutions  of  the  KZK  nonlinear  parab<^c  equation,  which  is  subject  to  the  same 
general  restrictions  as  NPE.  The  test  case  involved  propagation  in  an  air-filled 
rectangular  duct  with  rigid  walls  (duct  width  of  6.6  cm  in  z  direction)  in  which  a 
primary  wave  at  frequency  4350  Hz  was  excited  in  the  lowest  two  modes.  The  two 
numerical  solutions  were  compared  with  each  other  and  with  experiment. 

Figure  8  shows  the  numerical  solution  of  Eq.  (17)  (solid  lines)  compared  with 
the  numerical  solution  of  the  KZK  equation  (dotted  lines)  and  experiment  (circles) 
for  the  (a)  primary  wave  and  (b)  nonlinearly  generated  second  harmonic  component 
due  to  a  sotirce  level  of  150  dB  (re  20  /iiPa).  Note  that  the  solution  of  the  KZK 
equation  reveals  phase  errors,  as  a  result  of  the  parabolic  approximation,  at  both 
n  ~  1  and  n  =  2,  although  the  overall  amplitudes  are  in  good  agreement  with 
Eq.  (17).  Note  the  good  agreement  of  Eq.  (17)  with  the  measurements. 

Figiue  9  shows  comparisons  of  waveforms  due  to  a  source  level  of  154  dB.  The 
high  frequency  oscillations  in  the  first  and  third  coltunns  are  Gibbs  oscillations  due 
to  truncation  of  the  Foturier  series  at  20  harmonics.  Note  the  good  agreement  at 
all  four  ranges  between  Eq.  (17)  (first  column)  and  experiment  (second  column). 
However,  phase  errors  in  the  solution  of  the  KZK  equation  manifest  themselves 
as  differences  between  the  predicted  and  measured  waveforms  at  z  s  2.11  m  and 
z  =  5.36  m. 


Figure  8:  Propagation  curves  for  (a)  the  primary  wave  and  (b)  the  nonlinearly 
generated  second  harmonic  component  in  a  waveguide.  Solid  curves  are  obtained 
from  Eq.  (17),  dotted  curves  from  equations  based  on  the  parabolic  s^proximation, 
and  circles  are  measiirements.J 
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Figure  9:  Calculated  waveforms  based  ou  Eq.  (17)  (first  column)  and  on  the  para¬ 
bolic  wave  equation  (third  column).  Measured  waveforms  are  in  the  second  column. 
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